arXiv:chao-dyn/9706001vl 30 May 1997 


Quadratic Volume Preserving Maps 


Hector E. Lomelf 
James D. Meissfi 

Department of Applied Mathematics 
University of Colorado 
Boulder, CO 80309 

May 16, 1997 


Abstract 

We study quadratic, volume preserving diffeomorphisms whose inverse is also quadratic. 
Such maps generalize the Henon area preserving map and the family of symplectic 
quadratic maps studied by Moser. In particular, we investigate a family of quadratic 
volume preserving maps in three space for which we find a normal form and study 
invariant sets. We also give an alternative proof of a theorem by Moser classifying 
quadratic symplectic maps. 

AMS classification scheme numbers: 

34C20,34C35,34C37,58F05,70H99 

1 Introduction 

The study of the dynamics of polynomial mappings has a long history both in applied and 
pure dynamics. For example, such mappings provide simple models of the motion of charged 
particles through magnetic lenses and are often used in the study of particle accelerators 
[1]. Moreover, the quadratic, area preserving map, introduced by Henon [|j, is one of the 
simplest models of chaotic dynamics. 

Henon’s study can be generalized in several directions. For example, Moser ||] studied 
the class of quadratic, symplectic maps, obtaining a useful decomposition and normal form. 
Here we do the same for more general class of quadratic, volume preserving maps, with one 
caveat as we discuss below. 

Just as symplectic maps arise as Poincare maps of Hamiltonian flows, volume preserving 
maps are obtained from incompressible flows, and as such have application to fluid and 
magnetic field line dynamics HI- Moreover, one can argue that computational algorithms 
for flows should obey the “principle of qualitative consistency” j|]: if a flow has some 
qualitative property then the algorithm should as well. For the case of Hamiltonian flows 
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this leads to the construction of symplectic algorithms. A volume preserving algorithm 
should be used for a volume preserving flow, such as the motion of passive particle in an 
incompressible fluid || ||, ||, [l(| . 

Some of the properties of symplectic maps generalize to the volume preserving case. 
For example, a map that is sufficiently close to integrable and nondegenerate in a certain 
sense has lots of codimension one invariant tori 0, |T2| . PC , which are absolute barriers 
to transport [14]. Also, a perturbation of a volume preserving map with a heteroclinic 
connection can have an exponentially small transversal crossing [15]. Finally, the Birkhoff 
normal form analysis can be used to study the motion in the neighborhood of fixed points 

0 

Another motivation for the study of volume preserving maps is that they can be used 
as simple models for the study of transport in higher dimensions. The general theory of 
transport is based on a partition of phase space into regions between which transport is 


restricted by partial barriers [17j. For example, in two dimensions a partition is formed from 
intersecting stable and unstable manifolds of a saddle periodic orbit. In higher dimensions 
an analogous construction requires the existence of codimension one manifolds that separate 
the space [|l^]. In most cases it is difficult to find a dynamically natural construction of 
such manifolds; however, such manifolds do appear in volume preserving maps, and this 
leads easily to the construction of partial barriers. 

The computation and effective visualization of invariant manifolds in higher dimensional 
maps is itself an interesting problem [jl9[. In this paper we will study the intersections of 
the two dimensional stable and unstable manifolds in M 3 . 

Polynomial maps are also of interest from a mathematical perspective. Much work has 


been done on the “Cremona maps,” that is polynomial maps with constant Jacobians [20]. 
An interesting mathematical problem concerning such maps is the conjecture proposed by 
O.T. Keller in 1939: 


Conjecture 1.1 (Real Jacobian Conjecture) Let f 

Then f is bijective and has a polynomial inverse. 


be a Cremona map. 


This conjecture is still open. It is known that injective polynomial maps are automati¬ 
cally surjective and have polynomial inverses []2l], so it would suffice to prove that / is 
injective. It is easy to see (cf. lemma B below) that for the quadratic case, the condition 
of volume preservation implies injectivity, thus the Jacobian conjecture holds for quadratic 
maps. 

Even if the conjecture is true, the degree of the inverse of a Cremona map could be 
large. For example, the upper bound for the degree of the inverse of a quadratic map on R n 
is known to be 2" -1 (2^[. Thus in two dimensions the inverse of a quadratic area preserving 
mapping is quadratic, as was discussed by Henon [23, 201. More generally, Moser showed 
that quadratic symplectic mappings in any dimension have quadratic inverses ||. 

Henon found the normal form for the quadratic Cremona mapping in the plane. In 
this paper, we will correspondingly find the normal form for higher dimensional cases, but 
we assume that the quadratic, volume preserving mapping has a quadratic inverse (it is a 
“quadratic automorphism”). We give a complete classification of these diffeomorphisms. 

Such maps can be written as the composition of an affine volume preserving map and 
a “quadratic shear.” We give necessary and sufficient conditions for such shears to have 
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a quadratic inverse. As a first application of this concept, we give a simple proof of the 
theorem of Moser Q for the symplectic case. 

We also show that the quadratic automorphism in M 3 can be reduced to one of three 
normal forms. The generic case has four parameters: two govern the linearization of the 
map-the trace and second trace of the Jacobian matrix. The remaining two parameters 
determine the nonlinear terms represented by a single quadratic form in two variables. 

2 Quadratic Shears 

In this section we will study maps of the form 

x i—>• x + ^Q(x) 

where Q is a vector of quadratic polynomials. Throughout this paper we will write vectors 
of quadratic polynomials using the form Q(x) = M{x)x where M : M" —* GL(R") is a linear 
function into the set of n x n matrices that satisfies the symmetry property M (. x)y = M(y)x 
so that D x (M(x)x) = 2 M{x). 

Definition 2.1 We say that f : M™ —> M” is a quadratic map in standard form if f can be 
written as ^ 

f(x) = x + —M(x)x 

where M is a matrix valued linear function that satisfies M(x)y = M(y)x. 

It is important to notice that Df{x ) = I + M{x). 

Lemma 2.1 Let f(x) = x + |M( x)x be a quadratic map o/l n in standard form. The 
following statements are equivalent 

i) For all x € M n , det (Df(x)) = 1. 

ii) f is bijective with polynomial inverse. 

in) [M(x)] n = 0. 

Proof: 

We will show iii)=Mi)=M)=>iii). 

iii)=Mi) The condition implies that the matrix / + M (x) is invertible with inverse 
I — M(x) + M(x) 2 — ■ ■ ■ — (—1 ) n M{x) n ~ 1 . We notice that we can write 

f{x) - f(y) = (l + M (* - v)- 


So the function is injective. Using theorem A in [21], we conclude that / is bijective with 
a polynomial inverse. 




2 QUADRATIC SHEARS 


4 


ii)=^i) In principle, det(Df(x)) and det(Df 1 (/(x))) are polynomials in xi, x 2 , ■ ■ ■ ,x n . 
However, differentiation of f~ 1 (f(x)) = x gives 

det(-D/ _1 (/(x))) det(Df(x)) = 1, 

and therefore, since both are polynomials, det (Df(x)) has to be a constant independent of 
x. We notice that det (Df(x)) = det(D/(0)) = det(/) = 1. 

i)=>iii) Since det(/ + M(x )) = 1 and M is linear in x, then for any ( ^ 0 

det (M{x)-(I) = (-l) n Cdet {I + M(~x)) = (~1)”C"- 

This implies that the characteristic polynomial of M{x) is (—C)” an d therefore [M(x)] n = 0. 

■ 

At this point, we restrict to the case of quadratic maps in standard form whose inverse 
is also quadratic. We will see that the dynamics of such maps is essentially integrable, and 
is similar to the dynamics of a shear. We first establish the following characterization. 

Lemma 2.2 Let f(x) = x+^M(x)x be a bijective quadratic map o/M". Then the following 
statements are equivalent. 

i) f~ x is a quadratic map. 

ii) M(x) 2 x = 0. 

Hi) M(x)M(y)z + M(y)M(z)x + M(z)M(x)y = 0. 

k 

iv) f k (x) = x + —M(x)x for all k G Z. 

Proof: 

We will show i)=Mi)=Mii)=Mv)=M). 

i) =Mi) Let 

/ _1 (x) = x + ^N(x)x 

where N(x) is a matrix valued linear function that satisfies N(x)y = N(y)x. Then we have 
that 

111 1 

x = /(/” 1 (x)) = x + -N(x) x + -M(x + -N(x)x)(x + -N(x)x) 

= x + ^-N(x)x + l-M (x)x + \-M (x)N ( x)x + -M (N (x)x)N (x)x 
2 2 2 8 

This implies that N(x)x = —M(x)x and M(x)M{x)x = 0. 

ii) =^iii) By linearity and symmetry of M. 
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iii)=Mv) Let (Jk{x) = x + —M(x)x. Then 

l k l l 

gk{gi{x )) = x + -M(x)x + -M(x + -M{x)x){x + -M(x)x) 

1 i u ui up 

= x H- M(x)x -\ - M(x)M(x)x H- M{M(x))M{x)x 

2 2 8 

l + k kl 2 

= x H--— M(x)x -— [M(x)M(x)M(x)x + M(x)M(M(x))x\ 

2 8 

l T k , . 

= x H--— M(x)x. 

Therefore gk° gi = ■ On the other hand gi = f and < 7 0 = This implies that 


f k - 


iv)=H) Clear. ■ 

Definition 2.2 Let f : M" — » M n be given, in standard form, by f{x) = x + | M{x)x . If f 
satisfies any of the conditions of lemma [2A.J, we will say that f is a quadratic shear. 

A simple family of quadratic shears is determined by any vector v £ M" and a symmetric 
matrix P such that Pv = 0 according to M(x)y = ( x T Py)v , for then 

M(x) 2 x = (x T Pv)(x T Px)v = 0. 


We will see that, at least in the case n = 3, this is the most general quadratic shear. Moser’s 
normal form for symplectic, quadratic maps Q shows that the higher dimensional case is 
not quite this simple. From now on, we will concentrate on the special case n = 3. 


Theorem 2.1 A function f : R 3 —► M 3 is a quadratic shear in M 3 if and only if there is a 
vector v G l 3 and a 3x3 symmetric matrix P such that Pv = 0 and 


f(x) 


x + -(x T Px)v 


Proof: 

Since / is a bijection, we can define a new function g : S 2 —> S 2 . on the unit two 
dimensional sphere S 2 Cl 3 , in the following way. 


g{x) 


fix) 

!/(*)!' 


Using standard theorems of algebraic topology ^4|, we can argue that g has either a fixed 
point or an antipodal point (a point such that g(x) = —x). In any case, there is a constant 
K G M \ {0} and a vector x 0 ^ 0 such that f(x 0 ) = Kx 0 . We will show that K = 1. We 
notice that x 0 satisfies the following 


fix 0 ) = Kx o = x 0 + -M(x 0 )x 0 , 
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f 1 (Kx o) = x 0 = Kx 0 -— M(x 0 )x 0 , 

and therefore 

K 2 

K 2 (Kx 0 - x 0 ) = —M(x 0 )x 0 = Kx o - x 0 . 

We conclude that K satisfies K 3 —K 2 —K+l = 0 and hence I\ = 1 or K = — 1. Clearly, since 
/ is a bijection, f(^x 0 ) ^ 0 implies that M(x 0 )x 0 / —4x 0 . We conclude that f(x 0 ) ^ —x 0 
and therefore K = 1. 

It is clear that M(x 0 )x 0 = 0. Without loss of generality, we can assume that x 0 = e 3 = 
(1,0,0). Notice that M(e i) has to have the form 

M(e i) = ( 0 71 72 ) 

where 71,72 G M 3 . This fact, together with lemma |2.2(| , implies that the matrix M(e 1) 2 = 0, 
and therefore 71 and 72 are parallel. We can perform a linear change of coordinates and 
assume without loss of generality that / has the form 


/ Xi 

x 2 
\ x 3 


( X\ ^ 


^ a\ 

X 2 

+ X\X 2 

° 

V ) 




^ \ , 

— ^2 + ^2^3 

V I- 13 ) 

Let = M(ei),M 2 = M(e2) and M 3 = 
e 3 = (0, 0,1). It is easy to see that 



M(e 3 ) where e 3 = ( 1 , 0 , 0 ),e 2 


( 0 , 1 , 0 ) and 


Mi 


/ 0 a 0 \ 
0 0 0 
V 0 p 0 ) 


M 2 


/ 

a 

Hi 

Vi 


0 

H2 

V 2 

l 

P 

H3 

V 3 


and 



/ 

0 

v 1 

m 

m 3 = 


0 

Z ^2 

V2 


V 

0 

Z^3 

% 


To finish the proof, we need to show that the column vectors /it, v, rj and (a, 0, (3) of 
Mi, M 2 , M 3 are parallel to each other. We will show step by step that several of the entries 
are zero. We have two cases. 


(3 7 ^ 0. Using lemma [2.2] we conclude that 2M 2 e 3 + M 3 M 3 e 3 = 0. This implies that 
?y 2 = 0. We also have that M 3 = 0, so v 2 = 0 and r] 3 = 0. The condition M|e 2 = 0 
implies that fi 2 = 0 and this together with = 0 implies that v 3 = —a. Using the 
equation M 2 M 3 e 3 + 2M 2 e 2 = 0 we find that r/i = 0. Using M 2 = 0 and M 2 e 2 = 0, we 
find that the column vectors of M 2 are parallel, and the rest is clear. 
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• (3 = 0. The condition Af| = 0 implies that a = (1. v ?J = — y 2 and + v 2 y 3 = 0. 
The condition M| = 0 implies that ?y 3 = — v 2 and z/| — /U 2 7y 2 = 0. On the other hand, 
M|e 2 = 0 implies that /ii/r 2 + /i 3 ^i = 0 and M 2 e 3 = 0 implies that iq y 2 — rjiv 2 = 0. 

So it is enough to show that /qi ' 2 — n 1 y, 2 =0 and vpu 2 — ?/i/r 2 = 0. Clearly, if r/ 2 = (1 
then [i 2 = 0 and we would be done. So, we can assume that r\ 2 ^ 0. 

If 772 7 ^ 0 then t]i/Jj 2 = t\iv\lr\- 2 = v\v 2 . If v 2 = 0 then /r 2 = 0 and we would be done. 
Assume that v 2 ^ 0 and rj 2 ^ 0. This implies that /Ui^ 2 = //i/x 2 7 / 2 /^ 2 = — /u, 3 z^i 772 / 7^2 = 
Hlvi 772/^1 = ^ 2 ^ 1 - ■ 


3 Quadratic Symplectic Maps 

In this section we use the characterization of quadratic shears in lemma ||2.2H to give an 
alternate proof of the result of Moser ||] for quadratic symplectic maps. Recall that a map 
/ is symplectic if u)(Dfv, Dfv') = co(v, v') for all vectors v, v' G R 2 ” where lo is the standard 
symplectic form u>(v,v') = v T Jv' and J is the 2 n x 2 n matrix, 



Theorem 3.1 Let f be a quadratic symplectic map of (M 2n ,o;). Then f can be decomposed 
as f = ToS where T is affine symplectic and S is a symplectic quadratic shear. Furthermore, 
if S is any symplectic quadratic shear, then there is a symplectic linear map A such that 
A o S o A _1 (g,p) = (q + W(p),p). 

Proof: 

Let b = /(0) and L = Df( 0). Clearly L is a symplectic matrix and if we let T(x ) = Lx+b 
then S = T o / is a symplectic quadratic map in standard form, i.e. S(x) = x + \M(x)x , 
where M(x) is linear in x and satisfies the symmetry property M{x)y = M{y)x. Then S is 
symplectic providing 

(/ + M(x)) T J(I + M(x)) = J . 

Homogeneity of M(x ) implies that 

(1) M{x) T J = J T M(x) , 
and 

(2) M(x) T JM(x ) = 0 . 

Using (1) in (2) gives 0 = M(x) T JM(x) = J T M(x)M(x), and since J is nonsingular this 
implies 

(3) M{xf = 0 . 


Then lemma [2.2] implies that S’ is a quadratic shear. 



4 NORMAL FORM IN M 3 


To finish the proof, we follow Moser [[?] and define the null space of M(x) in the following 
way 

N = N(M) = {ye M 2n : M(x)y = 0,Va; G M 2n } = {y G M 2 ” : M(y) = 0} . 

Recall p5| that the u— orthogonal complement of a subspace E C R 2n is defined by E 1 - = 
{u G M 2ri : u(v,v') = 0,W G E}. We will show that N 1 - C N. For that purpose, we will 
use the following fact 

(4) M(z)M(x)y = M{x — y) 2 z = 0 

that follows from lemma |2.2(| , linearity and symmetry. 

Let u G N 1 - and x G M 2n . Now for any y G M 2n , (4) implies that M{x)y G N. 
Therefore u(y, M{x)u) = y T JM(x)n = —y T M(x) T Ju = —cu(M(x)y,u) = 0. This implies 
that M{x)u = 0 and hence u G N. Standard theorems in symplectic geometry (cf. [p5f ) 
imply that it is possible to find a Lagrangian space F such that N- 1 C F 1 - = F C N and a 
symplectic linear transformation A such that 


A (F) = {( q,p ) G R n X I" : p = 0}. 

Clearly, if S(x) = I + ^M(x)x is a symplectic quadratic shear, then so is 5 = A o S o A^ 1 . 
Assume that S(x) = I + \ M(x)x. Then A (F) C N{M). Therefore for all (q,p) Gl"x M n , 

M{q,p){q,p) = M(q,p)(0,p) = M(0,p)(q,p) = M(0,p)(0,p) . 

Since, in general, the matrix M(0,p) can be written in n x n blocks as 


M(0,p) 


A(p) B{p) ^ 
C(p) D{p) ) 


then M{0,p)(q,0) = 0 implies A(p) = C(p) = 0. Moreover, (1) implies D(p) = 0 and 
B(p) T = B{p). Thus finally we see that 

M(q,p)(q,p ) = (B(p)p, 0) 

where B{p)p is a gradient vector field. H 


4 Normal Form in M 3 

In this section we give normal forms for a quadratic diffeomorphism / of M 3 that preserves 
volume and has a quadratic inverse. Now lemma [f2.2| implies that if we let b = /(0) and 
L = Df{ 0), and T(x) = Lx + b, then the map S = T _1 o / is a quadratic shear. Then 
theorem 11 implies that S is of the form S(x) = x + |( x T Px)v where v G M 3 and P is a 
symmetric matrix such that Pv = 0. Depending on the relation between L and v, we have 
three cases possible cases; these can by distinguished by considering the space 


Z(v,L) = span{v , Lv , L 2 u} . 
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Theorem 4.1 Let f : R 3 —> R 3 be a quadratic volume preserving diffeomorphism. Then f 
can be written as the composition of an affine map T and a quadratic shear S, f = T o S, 
where S(x) = x + \(x T Px)v, v € M 3 and P is a symmetric matrix such that Pv = 0. 
Moreover , / is affinely conjugate to one of three possible normal forms, depending on the 
dimension of Z(v, L ) : 

i) dim Z(v,L) = 3. The map f is conjugate to 


(5) 


1 a + rx - cry + 2 + Q(x,y) \ 

x i 

V y ) 


where r and a are the trace and second trace of L, and Q(x,y) = ax 2 + bxy + cy 2 is 
a quadratic form. 


ii) dim Z(v,L) = 2. The map f is conjugate to 


( x 0 + ax + y+ Q(x, z) 
Vo ~ Px 
\ Zo + jz 


Hi) dim Z(v,L) = 1. The map f is conjugate to 


1 x 0 + ax + Q(y,z) 

y° - 

\ zo + y + fdz 


Proof: 

We know that / = L(x + \(x T Px)x) + b, and Pv = 0. To obtain the first normal 
form, perform a linear change of coordinates, x = Uf. Since the vectors v, Lv, and L 2 v are 
linearly independent, the transformation U can be defined by the following equations 

(5) = e 3 Ue 3 = v 

(6) U~ 1 Lv = ei Uei=Lv 

(7) U~ 1 L 2 v = e 2 + rei Ue 2 = L 2 v — TLv 

where, as we will see below, we will choose r = Tr(L). In the new coordinates the map 
becomes 


r =^-V(c/(e)) 

=C/- 1 6 + U-'LUf + i (f T U T PUf) U~ l Lv 
Ho + U-'LUtj + eM^if) 

where Q(^i 5 ^ 2 ) = \ PUf 2 ). Note that Q(f,ef) = ^ (f T U T Pv) = 0, so in the new 

coordinates the quadratic terms depend only on the first and second components. Moreover 
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in this coordinate system 

U~ l LUe\ =U~ 1 L 2 v = e 2 + re, 
U~ l LUe 2 =U~ l ( L 3 v - tL 2 v ) 
U~ 1 LUe 3 =U~ x Lv = e x 


The second equation can be simplified by noting that the characteristic equation for the 
matrix L is satisfied by L itself, and so L 3 — tL 2 + aL — 1 = 0 , where r = Tr(L ) and 
a = Tr 2 (L), the “second trace” of the matrix L, thus we get U~ 1 LUe 2 = t/ -1 (I — aL)v = 
e 3 — ae\. Thus we obtain 


TJ-'LU = 


Upon reverting to (x, y, z) as the names for the coordinates we get 




( x ° \ 


( tx - cry + z + Q(x,y ) \ 

u^nuix)) = 

o o 

+ 

x j 

l y ) 


To simplify this map further, we can conjugate, using the translation 


(,x, y, z ) ^ (x,y + y 0 , z + y 0 + z 0 ), 


to a map with x 0 = a, y 0 = 0 and z 0 = 0. This is the promised form. 

For the second case, assume that L 2 v = aLv — (3v , for some nonzero a and (3. This 
implies that the characteristic polynomial for L factors as (L—1//3I)(L 2 —aL+(31) = 0, and 
therefore, since L is nondegenerate, there exists a vector w ^ Z(v,L ) such that Lw = jjiv. 
We define the following change of coordinates. 


(8) 

U _1 u = e 2 

II 

(M 

03 

b 

(9) 

U~ 1 Lv = ei 

Uei = Lv 

(10) 

U~ 1 w = e 3 

Ue 3 = w 


As before, we note that in the new coordinates the quadratic term satisfies Q(e 2 ,£) = 0, so 
in the new coordinates the quadratic terms depend only on the first and third components. 
Moreover in this coordinate system we obtain 


U~ l LU 


fa 1 0 \ 
-(3 0 0 

V 0 0 i) 


This implies the form for the second case. 

For the third case, assume that Lv = av. We notice that there exist a vector w (j 
Z(v, L) such that Z(iv. L ) 0 Z(v, L) = M 3 . In fact, we can also find a constant (3 such that 
L 2 w — (3Liu + W = 0. We define the following change of coordinates. 


(11) 

U~ 1 v = c. 

& 

II 

03 

b 

(12) 

U~ 1 w = e 2 

Ue 2 = w 

(13) 

U~ 1 Lw = e 3 

Ue 3 = Lw 



5 DYNAMICS 


11 


As before, we note that in the new coordinates the quadratic term is Q(e i,£) = 0, so in 
the new coordinates the quadratic terms depend only on the second and third components. 
Moreover in this coordinate system we obtain 


U~ X LU 


(a 0 0 \ 

0 0 -i 

V 0 1 p ) 


This implies the form for the last case. H 


5 Dynamics 


The dynamics of the second and third cases of theorem 
ii), the z dynamics decouples from the (x,y) dynamics. 


[4.11 are essentially trivial. In case 
There are four special cases: 


• |/31 ^ 1. The plane z = (3z 0 /((3 — 1) is invariant, and is either a global attractor 
(|/3| > 1) or repellor (|/3| < 1). On the plane the dynamics is linear. 

• (3 = 1, z 0 ^ 0. All orbits are unbounded. 

• /? = 1, z 0 = 0. Every plane z = c is invariant, and the dynamics reduces to a two 
dimensional area preserving Henon map on each plane. 

• /? = — 1. Each plane z = c is fixed under f 2 . Restricted to a plane, f 2 is the 
composition of two orientation reversing Henon maps. 


For case iii) the ( y , z) dynamics is linear and decouples from the x dynamics. Generically, 
there is an invariant line on which the dynamics is affine. The invariant line can have any 
stability type. 


5.1 Generic case 


Equation (5) is the only nontrivial case. In general this map has six parameters, one from the 
shift, two from the linear matrix (the two coefficients of its characteristic polynomial) and 
the three coefficients of Q. However, generically, two of these parameters can be eliminated. 

Write the quadratic form as Q(x , y) = ax 2 + bxy + cy 2 . Generically a + b + c ^ 0 and 
we can we can apply a scaling transformation to set a + b + c = 1. Similarly if b + 2c 7 ^ 0 
the parameter a can be eliminated using the diagonal translation 


(x, y, z) ^ {x + 7 , y + 7 , 2 + 7 ), 7 = a/(b + 2 c) 


In this way, we obtain the final, generic form 


(15) 


^ x' \ 

/ 

a + tx + z + ax 2 + bxy + cy 2 \ 

y' = 


X 


V 

y ) 


a + b + c = 1 


There are four parameters in the system. Even ifa + 6 + c = 0 and/or b + 2c = 0, then 
other normalizations can be chosen to eliminate two of the parameters in (5). We will not 
study these special cases. 
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5.2 Periodic orbits 

Generically we can assume that a + b + c = 1 for the quadratic form in (5). The map (5) 
has at most two fixed points 

(16) x = y = z = x± = i r + a ± \J {t — a) 2 — 4 c^j 

born in a saddle node bifurcation at (t — a) 2 — 4a = 0. The characteristic polynomial of 
the linearized map at the fixed points is 

A 3 - t\ 2 + s\ - 1 = 0 

where the trace t and second trace s are 

t± = t T (2a + b)x± 
s± = a — (2c + b)x± 

We notice that 

t± — s± = ±yj (r — a) 2 — 4a 

t± — s T = ±(a — c)\J \t — a) 2 — 4a. 

The corresponding eigenvalues are illustrated in Fig |l[ It is easy to see (using the 
symmetric polynomials) that there are two lines in (t, s ) space where the stability changes: 
the saddle-node line t = s corresponds to an eigenvalue 1, and the period doubling line 
t + s = —2 corresponds to an eigenvalue —1. At the point t = s = —1 where they cross 
the eigenvalues are necessarily (—1, —1,1). Note also that when — 1 < f = s < 3 there is a 
pair of eigenvalues on the unit circle. There are two other curves of interest in the stability 
plane-these correspond to a double eigenvalue Ai = A 2 = r, or 

2r + l/r 2 = f r 2 + 2/r = s 

This gives the two curves shown in Fig |l|. One has a cusp at t = s = 3, where we have 
the triple root A = 1. The second crosses the saddle-node and period doubling lines at 
t = s = — 1. These are the two codimension two points. 

A fixed point with a one dimensional unstable manifold is called type A and one with 
a one dimensional stable manifold is called type B. The saddle-node and period doubling 
lines divide the plane into quadrants which alternate between type A and B. 

Having a pair of fixed points one of type A and one of type B, has interesting consequences 
for our map. For instance, the two dimensional manifolds serve as partial barriers to 
transport. Generically, they intersect along a one dimensional manifold. We have computed 
numerically some pairs of two dimensional stable and unstable manifolds. As an example, 
see figure ||. 

We have noticed that varying the parameter makes the one dimensional intersection 
bifurcate. Further investigation in this direction is the subject of future papers and more 
complete treatment will appear elsewhere. 
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The two fixed points (16) are born on the line t = s and move to opposite sides of this 
line for (r — cr) 2 > 4 a (x + is always on the right side). If \a — c\ is small, they are on the 
same side of the period doubling line, so that one is type A and the other type B ; however, 
when this parameter is large enough they can be on opposite sides, and therefore of the 
same type. This is determined by the sign of 

Sj- T f-j- 4-2 = 2 T t 1 (T T 2(n — c)x-}- 


When a = c, we have t± = so that the eigenvalues of the two fixed points are reciprocal 
(see section 5.3 for the explanation of this). 

Remember that, generically and without loss of generality, we can assume that a = 0. 
Therefore, we can plot stability diagrams for different values of r and a. The stability 
diagram in the (r, a) plane for the a = c case is shown in Fig |j. A more general case is 
shown in Fig [|. 

Periodic orbits can be studied by converting the map into a third order difference equa¬ 
tion. Let (x t ,y t , z t ),t = 0,1,... be a trajectory of the map (5), then the map can be written 


as 


x t +i = a + TXt- ux t - 1 + x t _2 + Q{x t ,x t - 1 ) 


Now it is clear that if this system of n quadratic equations has a finite number of solutions, 
then there are at most 2". We can rule out this degeneracy is most cases 


Lemma 5.1 Suppose a and c are not both zero, and let y± be the two (possibly complex) 
solutions of Q(/i, 1) = 0. Then if gJfgdfT ^ 1 for some integer 0 < k < n, the number of 
fixed points of /” for the map (5) is at most 2”. 


Proof: 

The “nonlinear alternative’" 


[p6| asserts that the number of complex solutions, counted 
with multiplicity, of a system of n polynomial equations in n variables is precisely the 
product of the degrees of the polynomials providing the system of equations obtained by 
discarding all terms but those of the highest degree in each equation has only the trivial 
solution. For our case, the resulting system is 


Q(x t ,x t - 1 ) = ax 2 t + bx t x t -i + cx 2 t -i = 0, t=,l...n, x n = x 0 


If any one of the x t = 0, then they all are zero, unless a = c = 0. Otherwise the general 
nonzero solution to this system is x t = p±x t -i, where ay 2 + bp, + c = 0. Setting x n = x 0 
requires pL k + p}_ = 1, where k + l = n. I 

For example, typically there are at most four fixed points of / 2 , giving a single period 
two orbit in addition to the two fixed points of /. However, there could be infinitely many 
period two orbits or none when /U +/ u_ = 1, giving a = c or n(_ = 1, giving b = ±(a + c). As 
an example, when a = c = 6/2 and a + r + 2 = 0, the line (x, 5 — x,x) (<5 — x, x, 6 — x) 
has period two where 5 is defined by 


a + (1 + cr)5 + ad 2 = 0 
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5.3 Reversibility 

A map is reversible if it is conjugate to its inverse by a diffeomorphism h that is an involution, 
thus 

ho f = f ~ 1 o h , h 2 = I . 

Some of the quadratic maps that we have considered have a reversor. Assume the generic 
case a + b + c / 0 (or equivalently Q{ 1,1) ^ 0, where Q is the quadratic form given in (5)). 
It is easy to see that, if a = c, then the map (5) has a reversor given by 


h{x,y,z ) 


^ Z + 7J 

y + v 

y X + 7] 


where rj = t — cr/(a + b + c). 

Note that when / is reversible and has fixed points, then the two fixed points have 
reciprocal eigenvalues-so if one is type A, the other is type B. Moreover, if the eigenvalues 
are complex, then the rotation rates have the same magnitudes at the two fixed points. 


Lemma 5.2 Let f be a quadratic map in normal form (5). Assume, generically, that the 
quadratic form satisfies Q( 1,1)^0 and (r — a) 2 ^ 4 a Q( 1,1) Then f is smoothly reversible 
if and only if Q(x, y ) = Q(y , x). 


Proof: 

Without loss of generality, we assume that Q( 1,1) = a + b + c = 1 and (r — a) 2 — 4a^0. 
Extend the map to C 3 . The imposed conditions imply that the map / has exactly two fixed 
points in C 3 . Suppose / is reversible and has a fixed points x±, then it is easy to see that 
h(x±) are also fixed points. In addition, Df(x± ) is conjugate to Df~ 1 (h(x±)). 

Since there are two fixed points, either h{x±) = x± or h(x±) = x T . Now the eigenvalues 
are invariant under a diffeomorphism, so in the first case the eigenvalues of Df(x + ) must be 
the same as those of H/“ 1 (x + ), and similarly for x_. This can only happen when t± = 
but this implies x + = x_, which can not happen by assumption. 

We conclude that h{x + ) = x_. This implies that t± = s T , which gives a = c, or 
Q{x,y) = Q(y,x). 

The other direction is proved by a simple computation, as described above. H 

Reversibility simplifies finding orbits of a map. Orbits that are invariant under h are 
called symmetric, and as is easy to see, they must have points on the fixed set Fix(h) = 
{x G M 3 : h{x) = x}. In our case this is the line x + z = — y = —r)/ 2, so a numerical 
algorithm for finding for symmetric orbits involves a one dimensional search. 

Similarly, if the stable manifold of one of the fixed points intersects Fix(h ), then the 
intersection point is on a heteroclinic orbit, for suppose x G Fix(h) 0 W s (x + ) , then x G 
W u (x + ), because hff n (x)) = f~ n (h(x)) = f~ n (x), so 

lim f n (x) = x + => lim h(f n (x )) = lim f~ n (x) = h(x + ) = x_ 

n —»oo n—>oo n—>oo 

Furthermore, suppose the stable manifold is two dimensional, and has the normal vector h 
at a point on Fix(h), then Dh{x)h is the normal to the unstable manifold at this point. 
This implies that the curve of heteroclinic orbits is tangent to the direction h x Dh{x)h. 



5 DYNAMICS 


15 


5.4 Bounded orbits 

For the Henon map, it is well known that the set of bounded orbits is contained in a square. 
For the volume preserving case, we will show that an analogue of this result also holds 
providing the quadratic form Q is positive definite: 

Theorem 5.1 If Q is positive definite then there is a re > 0 such that all bounded orbits 
are contained in the cube {( x,y,z ) : |x| < re, \y\ < re, \z\ < re}. Moreover, points outside the 
cube go to infinity along the +x axis as t —+ +oo or the -z axis as t —> —oo 

Proof: 

We start by writing the map in third difference form as 

x t+ i = a + Tx t - axt-i + x t -2 + Q{x t ,x t _ i) . 

Recall that a quadratic form Q(x, y) = ax 2 + bxy + cy 2 is positive definite iff a > 0, c > 0, 
and d = ac — b 2 /4 > 0. We will use the bounds obtained from completing the square: 

Q(x,y) = -x 2 + c(y + bx/2c) 2 > -x 2 , 

c c 

= (x + by/2a) 2 + -y 2 > -y 2 . 

a a 

There are three cases to consider, depending upon the relative sizes of x t , x t ~\, and .x t _ 2 : 

• l-Ttl > max(\x t -i\, |x t _ 2 |). The difference equation then gives 

x t +1 > Q{x t ,x t -i)-\a\-\Tx t \-\ax t _i\-\x t _ 2 \ 

> -xf - (\t\ + \a\ + l)\x t \ - \a\ , 
c 

Now since d/c > 0 there is a constant rei > 0, depending on a,r,a,a,b , and c such 
that when \x t \ > re l5 we have 

-x 2 - (|r| + |cr| + l)\x t \ - |a| > \x t \ , 
c 

In this case, we have x t +i > \x t \. Noting that we then have x t +i > \x t \ > \x t -i\, we 
can recursively apply this result to show that the sequence 


x t+k > x t+k -1 > . . . > \x t \ > Ki 


is monotone increasing. In fact, this sequence is unbounded; otherwise it would have 
a limit x t —> x* > rei, and this point would have to be a fixed point x = y = z = x* 
of the map. However, there are at most two such points, x±, and a simple calculation 
shows that rei > x±, so both fixed point points are excluded. 

• jxt— 2 1 > max(\x t \,\x t -i\). Inverting the difference equation and shifting t by one 
yields 


x t _ 3 =x t -a- Tx t -i + crxt -2 ~ Q(x t -i,x t - 2 ) ■ 
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Thus we have 

x t -3 < ~~Xt_ 2 + \ x t\ + \a\ + \rx t -i\ + \<7x t - 2 \ 

a 

< -~x 2 t _ 2 + (|r| + |cr| + l)|x t _ 2 | + \a\ 

a 

< ~\x t - 2 | , 

when |x t _ 2 | > « 2 , for a constant k 2 chosen as before, but with d/c replacing d/a. 
This implies that the sequence x t -k < x t _ fc+ i < ... < — |x t _ 2 | is monotone decreasing, 
negative, and unbounded. 

• \xt-i\ > max(\x t \, |x t _ 2 |). In this case we will see that the orbit is unbounded in both 
directions of time. For the forward direction, note that 

x t+i > ~Xt_! - (|r| + |cr| + l)|x t _i| - |a| 

a 

> \x t -\\ , 

when |x t _i| > k 2 . Thus x t+ \ > |x t _i| > \x t \, which is the situation covered by (i), and 
we get a monotone increasing sequence, providing x t +\ > Alternatively, note that 

x t -3 < -~x^_ 1 + (\r\ + \a\ + l)\x t -i\ + \a\ 

c 

< -\xt-i\ , 

when |x t _i| > k x . This gives x t _ 3 < — |x t _i| < — \x t - 2 \, so we are in the situation 
covered by (ii), which implies that the sequence approaches —00 providing |x t _ 3 | > n 2 


In conclusion, we have shown that an orbit is unbounded either as t —> ±00 providing it 
contains a point x t such that \x t \ > max{n 1 , n 2 ) = k. Note that n-y is a monotone decreasing 
function of d/c, so therefore we can define k by using d/max(a, c ): 


max(a, c) 


K = 


max(a, c ) 


2d 


|t| + |cr| + 2 + \ (|t| + |cr| + 2) 2 + 4 


a 


Finally, we investigate the asymptotic direction of an unbounded orbit. Recall that 
y t = x t -\ and z t = x t - 2 . Suppose that \x t \ > \y t \ > \z t \ > k, then each of the variables is 
eventually positive, so the orbit moves to infinity in the positive octant. Moreover, once all 
components are positive, we have 

Xt +1 _ QjxuVt ) a + Tx t — ay t + z t 

/y* Of* Of* 

Xt Xt Xt 

d led 

> -X t - ( T + cr + 1)-OO . 

c x t 

So the ratios y t /x t = x t _\/x t and z t /y t = x t _ 2 /x t _i go to zero, and the orbit approaches 
the positive x-axis as t —> 00 . Similarly if \z t \ > \yt\ > \x t \ > k, then eventually all 
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the components are negative, and so the orbit moves to infinity in the negative octant as 
t —> —oo. Once all components are negative, we have 


Zj-l 

z t 


< 


x t -3 
X t -2 

Q{yt,zt) x t — a — Ty t + az t 

z t Zt 

d I, .. . . . IckI 

- \z t \ + ( T + (J + 1) + 7 —r —> -OO . 

a \z t \ 


This implies the orbit moves to oo along the negative z axis. I 


6 Conclusions 


We have studied a family of volume preserving maps with the property that all entries 
are quadratic polynomials. We showed that these conditions imply that such maps are 
polynomial diffeomorphisms. Then we restricted ourselves to quadratic maps whose inverse 
is also quadratic. The class of maps studied is related to an old conjecture about polynomial 
maps called the Jacobian conjecture. 

A definition of quadratic shears was introduced and a characterization was given in 
general. A further characterization in three space was applied to find a normal form for the 
family. In three-space, the form of the generic case is similar in form to the area preserving 
Henon map and, generically, the map has two fixed points that can be either type A or type 
B. 


In addition, using our definition of quadratic shear and its characterization, we were 
able to give a simpler proof of a theorem of Moser classifying quadratic symplectic maps. 

Once we have found the normal form for the quadratic volume preserving map, there 
remain many enticing open problems. For example, we plan further computations to visu¬ 
alize the stable and unstable manifolds of the fixed points. Often these manifolds intersect, 
enclosing a ball; however, this is not guaranteed. Moreover, the heteroclinic intersections, 
which are generically curves, can fall in many homotopically distinct classes. We suspect 
that there are bifurcations between these classes, and that which occurs will depend, for 
example, on the complex phase of the eigenvalue of the associated fixed point. Heteroclinic 
orbits can be found most easily for the reversible case, as an intersections should occur on 
the fixed set of the reversor 

Another problem of interest is to obtain a characterization of quadratic shears in higher 
dimensions similar to the one we obtained in three dimensions. At this point, normal forms 
could be obtained using techniques similar to the current paper. 

Finally, as we discussed in the introduction, one of our main motivations for characteriz¬ 
ing the quadratic volume preserving maps is to study transport. If the two fixed points have 
disparate types, and their two dimensional manifolds intersect on a circle, then transport 
can be localized to “lobes” similar to the two dimensional case [^]. However, as Fig. |] 
shows, the intersections can be curves that spiral from one fixed point to the other. We 
plan to characterize transport for such cases. The existence, for the definite case, of a cube 
containing the bounded orbits (cf. theorem 5.1) will prove useful in this study. 
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7 Figures 



Figure 1: General stability diagram for a volume preserving map. 
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Figure 2: Two dimensional stable and unstable manifolds for the parameters 

a = c = 0.5, b = 0.0, a = 0.0, r = —0.3 and a = 0.0. 
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Figure 3: Stability Diagram for the reversible case a + b + c = 1, a = c and a = 0. There 
are no fixed points in the shaded region. The complex eigenvalues for the fixed point x + 
are shown. Those of x_ are reciprocal to those of x + . 
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Figure 4: Stability Diagram for a = —1/2,6 = l,c = 1/2, a = 0. Solid lines correspond to 
changes in the stability of x+, while dashed lines to x_. 













